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1 Field Theory 



These lecture notes want to illustrate the close connection between statistical mechanics 
and field theory not only on the formal level, i.e. that many concepts of one area can 
easily be taken over to the other one, but also on the level of actual calculations. To this 
purpose, the last section will demonstrate that a special statistical system, the binary 
fluid system, can be described by field theory in its critical behaviour. 

1.1 Canonical Formalism 

Let us briefly review the setup of the canonical formalism of quantum field theory. For 
simplicity, we restrict ourselves to the case of a Hermitean scalar field <j>{x) in Minkowski 
space. 

There is a Hilbert space TC of physical states containing the vacuum |0) as the 
state of lowest energy. We assume that every state vector can be written as a linear 
combination of products of field operators acting on the vacuum |0). Furthermore, TC 
carries a unitary representation U(a, A) of the Poincare group, where a denotes a space- 
time translation and A = (A^^) a (homogeneous) Lorentz transformation. The vacuum 
is singled out by the fact that it is the only state vector which is invariant under U(a, A). 
The generators P M of space-time translations, 

U{a,l) =exp(iPX) > (!) 
have a spectrum which is confined to the closed forward light conef] ( "spectral condition") 

p 2 = Pfl p^>0; p°>0. (2) 

P° here denotes the Hamiltonian H. The scalar field <f)(x) is a Hermitean, operator-valued 
distribution on TC, which transforms covariantly under Poincare transformations: 

U(a,A)(f)(x)U j (a,A) =(j)(Ax + a) . (3) 

The theory is quantised by imposing canonical commutation relations between the fun- 
damental field 4>(x) and its canonical conjugate momentum ir(x) at equal times: 

[0(t,x),7r(t,y)] = i£( 3 )(x-y), 

[0(t,x),0(t,y)] = [7r(t,x),7r(t,y)] = 0. 1 ] 

Causality is guaranteed by requiring locality for the field, i.e. for space-like separations of 
the arguments the field operators must commute with each other: 

[<j>(x),<j>(y)] = for (x-y) 2 <0. (5) 

As the appearance of the 5-distribution in (^) indicates, local field operators as <p(x) 
must be viewed as operator-valued distributions rather than functions and should be 
smeared out with appropriate test functions, 

</>(/) = |d 4 x/(x)0(x). (6) 
1 We use the Minkowski metric = diag (1, —1, —1, —1) 
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The smearing of 4>(x) also in the time variable is necessary for interacting fields. Then, 
however, one can no longer postulate well-defined commutation relations at equal times 
An alternative to the Hilbert space formulation is to avoid operators and states and 
to introduce the functional integral as an independent way of quantising a theory. To 
make contact between both approaches, we introduce Green's functions, which play a 
fundamental role in both formulations. 

1.2 Green's Functions 

Vacuum expectation values of products of field operators 
W{x 1} ...,x n ) := (O|0(x 1 )...0(x n )|O) 

(7) 

= ( O|0(O, x^e- 1 ^-* 2 ^ . . . e- i (*"- 1 -*") H 0(O, x n ) |0 ) 

are called n-point correlation functions or Wightman functions. A more prominent role 
in the Minkowski space formulation of field theory is played by time- ordered Green's 
functions which are defined as the vacuum expectation values of time-ordered products of 
field operators, 

G(x h ...,x n ) := (O|T0(t 1 ,x 1 )...0(t n ,x n )|O) 

r i ( 8 ) 

= (0\T< 0(0, x 1 )e- i (* 1 -* 3 ) H . . . e- i (*"- 1 - i ^ H 0(O, x n ) Uo } ; 

here T denotes the Dyson time ordering symbol, 

T{0(xi) . . . 0(x n )} := 0(^(1)) . . . (j){x^ n )) with t n{1) > t n{2 ) > > t n{n) . (9) 

They can be directly related to S-matrix elements via the well-known Lehmann-Symanzik- 
Zimmermann (LSZ) reduction formulae. E.g., the S-matrix element between an out-state 
with momenta q\, . . . , qi and an in-state with momenta p%, . . . ,p n is given by 



>i,... ,p n ,in|S|gi, . . . , ft, out) = 

= (iZ- 1 / 2 )^ 1 [d% . . . d 4 x t exp ( i £ p k ■ y k - i E q r ■ xA 

J \ k=l r=l J 

x (D yi + m 2 ) ■ ■ ■ (a xi + m 2 ) G(yi, ..,,!,), 



(10) 



where Z is an appropriate normalisation factor. For more details on the LSZ formalism 
the reader is referred to a standard textbook on quantum field theory, e.g. Refs. @, [f|. 

The most important fact about Green's functions is that they contain all the informa- 
tion about the theory. So, given either all the Wightman or all the time-ordered correlation 
functions we may construct the Hilbert space and the unitary representation U(a, A) out 
of them (this is the essence of the so-called reconstruction theorem, @], fj). Thus, we may 
use the Green's functions to define the theory, thereby avoiding the operator formulation 
altogether. 

1.3 Euclidean Time Formulation 

Another important ingredient in our formulation of Euclidean field theory is the analytic 
continuation to imaginary times: 

x° -> -ix 4 , x 4 e R. (11) 
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This is convenient for several reasons. Generally speaking, it improves the analytic be- 
haviour of the various relevant functions. In perturbation theory it simplifies the calcu- 
lation of Feynman diagrams, e.g. because of the positivity of the energy denominators. 
Moreover, since in Euclidean space all distinct world-points are space-like separated, the 
Euclidean Green's functions are automatically symmetric (for bosons) in all their space- 
time arguments, and we do not need any time-ordering. Thus, in Euclidean space we have 
to deal with only one kind of Green's functions and both, the Wightman and the time- 
ordered Minkowski space Green's functions, may be obtained from the Euclidean ones by 
appropriate analytic continuation. Last not least, the functional-integral formulation to 
be introduced in section ( |1.5|) is much better defined in the Euclidean formulation due to 
the positivity of the Euclidean Green's functions and may be interpreted as a stochastic 
process. 

Consider the analytic continuation of the Wightman functions to complex time argu- 
ments x° — > x° — ixf , xf G 1R : 

W(x? - ixf, xi; . . . ; x° n - ia£, x n ) = 

(12) 

= ( O|0(O, Xl )e- i( - x i- x °^ H e-^- x ^ H (f)(0, x 2 ) . . . 0(0, x n ) |0 ) . 

Due to the positivity of the Hamiltonian this is well-defined and analytic at least for 

x\ > x\ > . . . > x* , (13) 

since we have exponential suppression in this domain. We define the Euclidean Green's 
functions, also called Schwinger functions, as the x® — > limit of the analytically continued 
Wightman functions: 

<S(xi, x±; . . . ; x n , x n ) := W(— ix^, Xi; . . . ; — ix n , x n ) = 

(14) 

= ( O|0(O, Xl ) e -(4-4)H . . . e-^-i-^^O, x n ) |0 ) . 

For the time being this is only legitimate if the condition (|13|) is fulfilled. It can, however, 
be shown [Q that the region of analyticity is much larger and the Schwinger functions are 
actually well-defined for all non- coinciding Euclidean world-points: 

Xi 7^ Xj for i 7^ j . (15) 

Let us demonstrate this for the two-point function: The translational invariance of the 
vacuum implies that W depends only on the difference of its arguments 

W(x l7 x 2 ) = W(x 1 -x 2 ), (16) 

which can be seen explicitly from (|i~2"D. As a function of complex x°, W is analytic in the 
lower half plane. For permuted arguments we have 

W(x 2 ,x 1 ) = W 7T (x 1 -x 2 ) = W(x 2 -x 1 ) with W n {x) :=W{-x) , (17) 

where W n is now analytic in the upper half plane according to (|TB|). For space-like 
separations, (xi —x 2 ) 2 < 0, however, the Wightman function is symmetric due to locality: 

W(x 1 ,x 2 )=W(x 2 ,x 1 ), (18) 
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Figure 1: The analyticity region of the two-point function W(x) in the complex 
x°-plane. The Wightman function W is obtained by approaching the real line 
from below. The Schwinger function is defined on the imaginary x°-axis. The 
time-ordered Green's function G is obtained from the Schwinger function by 
means of a Wick rotation, indicated by curved arrows. 

implying 

W v {x) = W(x) (19) 

in the non-discrete region 

- |x| <x° < |x| (20) 

of the real axis. The well-known edge-of-the-wedge theorem then guarantees that W and 
W n form a single, analytic function in the union of their domains, i.e. everywhere in the 
complex x° plane except along the cuts |x°| > |x|. 

The Schwinger functions and their properties have been studied in an axiomatic setting 
by Osterwalder and Schrader ||. For our purposes the most interesting properties are 

1. Euclidean covariance: 

S(Ax 1 + a,...,Ax n + a) =S(x ll ...,x n ) , A G SO(4) , (21) 
where the Lorentz group of ([|) is replaced by the group of Euclidean rotations. 

2. Symmetry: 

S{x a (i), X a (n)) = S(X!, ...,x n ), (22) 

where a is a permutation. This rests on the fact that distinct Euclidean points are 
always space-like relative to each other. 
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3. Reflection positivity: This property replaces the Hilbert space positivity and the 
spectral condition (Q) of the Minkowskian formulation and is necessary to guarantee 
that the Euclidean correlation functions may be continued back to Minkowski space. 
Formally, it is defined as 

/ d4 ^i • • • d 4 ^ d% . . . d 4 yj f*(x 1 , ...,x t ) fj(y u Vj) 
i,j J (23) 

xS(9x 1 ,...,9x i ,yi,...,y j ) > 0, 

where 9 is the Euclidean time reflection, 

fl(x,* 4 ):=(x,-x 4 ), (24) 

which, roughly speaking, replaces the Hermitean conjugation in Minkowski space. 

A Euclidean quantum field theory may be defined by the set of all its Schwinger 
functions. From them the Wightman functions in Minkowski space can be reconstructed 
as boundary values: 

W(xi, • • • ,x n ) = lim 5(xi,e a + ix?; . . . ;x n ,e n + ix° n ) , (25) 

£i->0, 
Si-e i+ i>0 

i.e. approaching the real x°-axis from below in the case of the two-point function, see 
Fig. |l[ The advantage of the Euclidean formulation, however, is that the Schwinger 
functions obey simpler properties and are easier to handle than Wightman functions or 
field operators. In particular their symmetry is the crucial property which opens the way 
to a representation in terms of functional integrals. 

The time-ordered Green's functions G(x\, . . . , x n ) may also be obtained from the 
Schwinger functions, namely by approaching the real a;°-axis through a counter-clockwise 
rotation of 7r/2 as indicated in Fig. IB, 

G( Xl , ...,x n )= lim 5(e iQ a;?, . . . , e ia x° n ) . (26) 

This rule, which obviously is simpler than the one for the Wightman functions, is called 
Wick rotation. 



1.4 Example: Free Field 

To illustrate the analytic continuation discussed above, let us consider the simplest case 
of a free massive Hermitean scalar field <f)(x), decomposed in momentum space as 0] 

m = I (44^ 

with k° = ujk = Vk 2 + m 2 . The annihilation and creation operators satisfy the commu- 
tation relations 

[a(k),a)(k')} = (27r) 3 2cj fc 5 3 (k-k'). (28) 



a(k)e-' lkx + a\k)e 



ikx 



(27) 
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The two-point Wightman function then is given by 



W{x)=W(x,0) = (O|#r)0(O)|O> = A + (x) = J-0- 4 5{k 2 - to 2 ) 9(k°)e- ik \ (29) 
The time-ordered two-point Green's function, the Feynman propagator, reads 

iA P (x) := 0) = ^(x°)^(x) + 9(-x°)W(-x) = jj-^ fc2 ^ - . £ ■ (30) 
Performing the analytic continuation of W(x) yields 

G(x, 0) is obtained from S(x, 0) by means of a Wick rotation of x 4 = ix° in a counter- 
clockwise sense. In the integral we have to rotate at the same time k 4 = ik° in a clockwise 
sense. The resulting integration path in the complex /c°-plane bypasses the poles at 
±Vk 2 + m? in the way indicated by the +ie prescription. It goes from +oo to — oo, and 
an additional sign change has to be introduced to bring it in the form above. Note that 
the long distance behaviour of the Schwinger function is exponential and governed by the 
mass to: 

^■°)'^2(2£^ eXp( - m|x|) - (32) 
1.5 Functional Integral Formulation 

The symmetry property of the Schwinger functions means that the Euclidean fields com- 
mute. In this respect they behave like classical fields. This allows us to interpret the 
0(x)'s as classical random variables instead of quantum mechanical operators. In this 
interpretation the Schwinger functions are the n— point correlation functions or moments 
of an appropriate probability measure d/x[0]: 

S(xi,...,x n ) = . ..(f>(x n ) ) = J d/i[0] <j>(xi) . ..<j)(x n ) . (33) 

We may define a generating functional for the Schwinger functions as the Laplace trans- 
form of the probability measure, 

Z[j] := ^exp d 4 x(p(x)j(x)^ ^ 

(34) 

00 1 r 

= Yl ~1 <^ x i ■ ■ -d 4 x n j(xi) . . .j(x n ) ( 0(xi) . . . <f>(x n ) ) , 

n=0 U - J 



from which the moments are obtained by functional differentiation with respect to the 
external sources j(x): 



Thus the generating functional Z[j] summarises all the information about the theory in 
a highly condensed form. 

Formally, one may decompose d/x[0] in a product of a "Lebesgue" measure T>[4>] on 
function space and a normalised weight factor 

d/i[0] = lx?[0]e- SE M, (36) 

where Se[4>] is the Euclidean action of the field <j>, 

SbM = / d 4 x {- 8^{x) + ^V(x) + V (0(x))| . (37) 

This yields the functional integral formulation of the Euclidean correlation function and 
the generating functional: 

= \ fv[<f>] exp(-S E [0]+ jd 4 x0(x)j(x)) , (38) 
= iyb^^.-.^arOe-^, (39) 
where -Z is introduced to normalise the generating functional to Z[0] = 1, 




However, as it stands, (^) is a purely formal definition since a translational invariant 
Lebesgue measure like T>[<j)} does not exist on the infinite-dimensional function space of 
the continuum field theory. 



1.6 Lattice Regular isat ion 

We can give a precise meaning to the functional integral expressions of the preceeding 
section, if we refer to a lattice regularisation of the Euclidean field theory. This means, we 
replace the Euclidean space-time continuum by a finite hyper-cubical lattice with lattice 
spacing a, 

x^ = an fJi , n M integer- valued . (41) 

One may think of the lattice field <f)(x) as the average of the continuum field over the 
volume of the elementary cells of the lattice. In momentum space, the lattice regularisation 
reflects itself in the restriction of the momenta to the first Brillouin zone, 



K <P,<-\- (42) 



a a 



On a finite lattice the function space is finite-dimensional, and the measure T>[(p] can be 
defined as the multi-dimensional Lebesgue measure 

P[0]:=nd0(x), (43) 
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where x runs over all lattice points. Space-time integrals are replaced by sums and 
the differential operators <9 M are replaced on the lattice by nearest-neighbour forward 
differences, 

/d 4 *-£" ; 



a 

X 



1 ( v (44) 

<9 M 4>(x) -> A^0(x) := - \ 6(x + afi) - </>{x)\ , 

where fi is the unit vector in //-direction. This yields the lattice analogue of ([37]): 

SM = -E« 4 {(^0(x))(A,0(x)) + ^0 2 (x) + V(0(x)) 

r 4 / \ / m 2 n 2 \ n 2 ~\ ( 45 ) 

= 5> 2 {- E(*(*M* + «/*)) + (4 + ^J0 2 (*) + y V^))j . 

The lattice functional integral version of the Schwinger functions, 

( <j>{x x ) . . . <j>{x n ) ) = i I J] d0(x) 0(x x ) . . . 0(x n ) e- 5E ^ , (46) 



may also be derived from the expression ( |T4| ) by discretising the Hamiltonian, breaking 



up the exponential into a product of evolution operators between two neighbouring time- 
slices and inserting complete sets of field eigenstates at every time-slice - in the same 
fashion as the path integral is derived in ordinary quantum mechanics, see e.g. Refs. 0, 

I, I- 

The continuum theory is obtained back from the lattice version by taking two inde- 
pendent limits: 

1. Thermodynamic limit: taking the lattice size and - with it - the number of 
degrees of freedom to infinity; 

2. Continuum limit: making, loosely speaking, the lattice structure vanish, i.e. by 
taking the lattice spacing a to zero while keeping physical quantities like the mass 
gap at their actual value. We will see later that the continuum limit of a lattice 
field theory corresponds to a statistical lattice model approaching one of its critical 
points. 

1.7 Analogy to Statistical Mechanics 

The functional integral formulation of the generating functional of Euclidean lattice field 
theory, 

Z\j] = — f T\d<f)(x) e -SvMx)]+fd±x<p{x)j(x) ? ( 47 ) 

Z J X 

exhibits a close analogy to the statistical partition function of a system of spin variables 
(or magnetic moments, in general: a local order parameter) 4>(x) on a crystal, coupled via 
next-neighbour interactions (see (|45|) ) and additionally to an external source j(x) (e.g. 
magnetic field). The probability weight e~ SE ^ corresponds to the Boltzmann factor e _/3H . 
The vacuum expectation value of the field 

(48) 
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QUANTUM FIELD THEORY 


CLASSICAL STATISTICAL MECHANICS 


scalar field <f>(x) 

generating functional Z 

Euclidean action Se 

vacuum expectation value ( <p(x) ) 

Schwinger function S(xi,X2) 

1 

inverse mass — 
m 

Lagrangean formulation in d dimensions 


spin variable, local order parameter <f>(x) 

partition function Z 

Hamiltonian f3H 

mean magnetisation M 

spin-spin correlation function ( 0(xi)0(x2) ) 

correlation length a£ 

Hamiltonian formulation in d + 1 dimensions 



Table 1: Analogy between Euclidean quantum field theory and statistical mechanics 



corresponds to the mean magnetisation M per site of a ferromagnet, and the two-point 
function 

(ftxMfa)) (49) 

is equal to the spin-spin correlation function. The (dimensionless) correlation length £, 
which describes the spatial extent of fluctuations in a physical quantity about its average, 
governs the exponential decay of the correlation function in the long distance limit, 

( 0(^i)0(^ 2 ) } ~ e-l* 1 - 9 *!/* . (50) 

It is related to the mass gap (inverse Compton length) m = E\ — Eq of the field 
by 

am 

as can be seen from the expansion (|I4"D 

(0(a;)0(O)> = |(O|0(O)|n)i 2 e-^-^), (52) 

n>0 

where we inserted a complete set of energy eigenstates and, for simplicity, have chosen 
x = 0. In the long (Euclidean) time limit the two-point function consequently decays 
exponentially 

(0(x)0(O)> ~ e"^ 1 "^ 4 . (53) 

As is well-known in statistical mechanics, the magnetic susceptibility is related to the 
correlation function by 

X = £(0(x)0(O)> (54) 

X 

and therefore equals the propagator at zero momentum. 

The analogy between Euclidean quantum field theory on a lattice and statistical me- 
chanics - summarised in Tab. |T| - has turned out to be very useful. Many concepts and 
methods of statistical mechanics have been applied to field theory, and conversely, the 
field theoretic renormalisation group is an important tool in statistical mechanics. We 
will make extensive use of this analogy in the second part of this talk. 
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2 Critical Phenomena 



2.1 Critical Points 



We now turn to the question of a continuum limit. If a continuum limit with a finite 
physical mass exists, it means that by a suitable choice of the bare parameters we can 
approach a limit where a goes to zero while m remains finite. According to ([51]) the 
(dimensionless) correlation length has to diverge in that case. 

A point in the coupling constant space - the space whose coordinates are the param- 
eters or coupling constants of a theory -, where £ diverges and where the first derivatives 
of the relevant thermodynamic potential, say the Gibbs potential, exist, is called a crit- 
ical point or a second order phase transition. For most systems, the behaviour of many 
quantities near a critical point is governed by simple power laws, the so-called scaling 
behaviour. As the temperature T approaches its critical value T c , the correlation length 
and susceptibility, for example, diverge according to 



T -T r 



X 



T — T r 



(55) 



with certain critical exponents v and 7. The relation ~ is to be understood in the sense 
that 

* lim^M = <x. 
The magnetisation in the low temperature phase vanishes like 



M 



T-T c 



P' 



(56) 



In four dimensions these laws are modified by logarithmic corrections. 



2.2 Universality 

The investigation of various systems of statistical mechanics near their critical points has 
revealed the property of universality. To be precise, the systems fall into a relatively small 
number of universality classes. The members of a class show identical critical behaviour 
in the sense that their critical exponents as well as certain other universal quantities 
are equal. As an example we mention that water at its triple point falls in the same 
universality class as a ferromagnet at the Curie point, which shows off in the same critical 
exponent v = 0.630. 

The universality classes are distinguished by 

1. the number of dimensions of space (or Euclidean space-time), 

2. the number of degrees of freedom of the microscopic field, 

3. the symmetries of the system. 

Universality means that the long-range properties of a critical system do not depend 
on the details of the microscopic interaction. In particular, also the size of the lattice 
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spacing becomes unimportant for the large-distance behaviour of correlation functions, if 
the correlation length is large. According to the scaling hypothesis, the correlation length 
is the only relevant length scale for the system near criticality. This hypothesis leads to 
various relations between critical exponents, such that only two independent exponents 
remain. 



2.3 Renormalisation Group 

Scaling theory and universality have found a theoretical basis in the Kadanoff- Wilson 
renormalisation group. Here we shall try to give only a brief sketch of the basic ideas 
of the Kadanoff- Wilson renormalisation group. For a more detailed presentation see e.g. 
Ref. 0. 

The original action S with a cut-off A (e.g. the boundary of the Brillouin zone, -) is 
considered to be embedded in an infinite-dimensional space of actions 

S = y £K i S ii (57) 

i 

where the coefficients Ki are called coupling constants and 

3 = £A (58) 

X 

with so-called local operators £j. They are functions depending on the fields at points x 
and a finite number of points near x, and having a certain engineering dimension, e.g. 
<f) n (x) has dimension n. A renormalisation group transformation R\ is a mapping 

R x : S i-> S {X) (59) 

in this space such that both S and describe the same physics at large distances, but 
the cut-off A gets lowered by a factor A > 1: 

A^-A. (60) 
A 

In other words, is obtained from S by integrating out degrees of freedom with high 
momenta near the cut-off. R\ can be described in terms of the changes of coefficients 

Ki^Kl**. (61) 

The most important points in the space of actions are the fixed points S*, 

R X S* = S* , (62) 

in particular those where the correlation length is infinite. Near a fixed point, the action 
of R\ can in general be linearised and diagonalised such that in a suitable basis 

K a = K* + 6K Q (63) 

it reads 

= K* a + X da 8K a . (64) 
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Figure 2: Vicinity of a fixed point S* in coupling constant space (shown are 
only two couplings K 1 , K 2 out of infinitely many) with one relevant (i.e. leaving 
the fixed point) and one irrelevant direction (approaching the fixed point). The 
action of a renormalisation group transformation R\ mapping an action S into 
another action is also shown. 



Those terms with negative scaling dimension d a get suppressed after repeated appli- 
cation of the renormalisation group transformation and are called irrelevant, since their 
presence does not affect the long-distance physics. The terms with positive d a , which are 
a few in general, are relevant. The values of the corresponding coefficients K a are decisive 
for the long-distance physics. Terms with d a = are called marginal. 

In this picture universality emerges in the following way. Two original actions S' and 
S", which belong to the domain of attraction of the same fixed point, are mapped under the 
action of the renormalisation group into the neighbourhood of the same low-dimensional 
manifold 

S = S*+ J2 K «S a , (65) 

relevant a 

where for simplicity we assume that no marginal operators are present. The critical 
behaviour is then determined only by the few relevant operators in the vicinity of the 
fixed point. In particular, it can be shown that the critical indices are simple algebraic 
combinations of the dimensions d a belonging to them. Thus the fixed points of the 
renormalisation group determine the universality classes of the actions. 

In four dimensions, perturbative calculations indicate that for the scalar field theory 
under consideration there are only two relevant operators, which are essentially the mass 
term <p 2 (x) and the linear term (j>{x), which appears when an external "magnetic" field is 
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present. The quartic self-interaction <j) 4 (x) is marginal, but its coupling g decreases under 
the renormalisation group transformation. The associated fixed point therefore has g = 
and belongs to a free field theory. It is called the Gaussian fixed point. Indications by 
non-perturbative methods lead to the result that the Gaussian fixed point is with high 
certainty the only fixed point for this theory. 



3 Interface Tensions of Binary Fluid Systems 

That universality and the renormalisation group are not just mathematical constructs 
but can be tested experimentally is demonstrated in the following example. It exhibits 
again the connection between statistical mechanics in d + 1 dimensions and field theory, 
especially the Euclidean functional integral formulation of it, in d dimensions. 



3.1 Phenomenology of Binary Fluid Systems 



Consider the mixing and separation of two fluids: Trying to mix Cyclo-hexane (CeH 12 ) 
and Aniline (C^H^N H 2 ), one notes that below a critical temperature of T c = 30.9°C, 
both fluids separate spontaneously into two pure phases consisting of Cyclo-hexane and 
Aniline respectively, one on top of the other. The surface tension r of the interface 
between them vanishes as the temperature T approaches T c , and one measures a scaling 
law for the reduced interface tension (k is Boltzmann's constant) a := r/kT, 



fj, = 1.26 ±0.01 



(66) 



t :-- 



T — T r 



o"o is the critical amplitude of the interface tension. "~" denotes the critical behaviour of 
the interface tension near T c , as defined earlier. 

Above T c , the two fluids mix perfectly, and hence a homogeneous phase is prepared. 
Approaching the critical temperature from above, the experimentalist notes that the 
mixed phase becomes less and less transparent, and at the critical point, the system is 
completely opaque, indicating that its correlation length £ + diverges (the " + " denotes 
approach of T c from above). No latent heat is set free, and the system hence exhibits a 
second order phase transition. The correlation length is measured to scale above T c like 



0.630 ±0.002 



(67) 



Below T c , the same critical behaviour is expected but with a different critical amplitude 
£o~ for the correlation length. (In this section £ is considered to be dimensionful.) 

There are many other binary fluid systems like Isobutyric acid and water, Triethy- 
lamine and water, and also systems of one fluid and one gas, which show the same be- 
haviour. Although the critical amplitudes o"o , £o vary considerably from system to 
system, the critical exponents agree within a few percent and obey Widom's scaling law 



fi = 2i/. (68) 

The scaling hypothesis indeed predicts /i and v to be universal and also gives Widom's 
scaling law. Renormalisation group calculations also predict the correct value for the crit- 
ical exponent v [D]. Binary fluid systems therefore seem to obey the scaling hypothesis, 
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and the critical exponents suggest that they lie in the same universality class as the three- 
dimensional Ising model. As a consequence of this, and although the critical amplitudes 
vary considerably from system to system, the dimensionless quantities 

R±:=<ro(to) 2 (69) 

should be universal, and experimentally this is indeed found to be the case: 

R e *v = 0.38 ± 0.02 . (70) 

Because £_ is hard to measure, R- is not easily accessible experimentally. In contradis- 
tinction, it turns out in field theory that the low temperature value is easier to obtain 
than R + . Monte Carlo data of Mon pf yield 

R^ c = 0.36 ±0.01 , (71) 

confirming the scaling hypothesis. In this calculation, finite-size effects on a have been 
shown to play an important role. On the other hand, a first field theoretic treatment gave 
an unacceptable value of R + = 0.20 [Q. Is there a real discrepancy between experiment 
and field theory? 



3.2 Field Theory of Binary Fluid Systems 

In the framework of field theory, one investigates critical phenomena of the systems dis- 
cussed above in the context of a Euclidean, massive and real $ 4 -theory, which is believed 
to be in the same universality class as the Ising model and the binary fluid systems. One 
may motivate this as follows: One needs a local order parameter which vanishes above 
T c and is nonzero below T c , indicating the strength of the symmetry breaking. The dif- 
ference between the concentrations of the two fluids A and B at a point is surely a good 
candidate, 

$(x)ocp A (x)-p B (x), (72) 

and symmetry breaking is - as in a ferromagnet - achieved spontaneously. 
One therefore considers in field theory the bare Lagrangean 



£ = i(9 M $)(^$) + V($) (73) 



with the double well potential 



Vm = -\ml& + + g = I S0 (* 2 - vtf , (74) 



whose minima in the phase with broken symmetry (vq > 0) lie at 

$ = ±V = ±< 



3mg 



\ 9o 



(75) 



One of the minima, $ = vq, is identified with the system being in the A-phase, the other 
one, $ = —Vq, with the -B-phase. In the phase with spontaneously broken symmetry, 
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a 

/ Kink 

" v o 



Figure 3: The kink solution in $ 4 -theory. 

T < T c , when the two components of the fluid separate, one therefore obtains the picture 
of Fig. |3| for $, when one moves on a trajectory perpendicular to the interface between 
the two fluids. The interface has been taken to be perpendicular to the ^-direction, the 
Euclidean time. It corresponds to the well-known kink in field theory, upon which will be 
dwelled in more detail below. 

Mind that mo , go are bare quantities and need renormalisation. The renormalised 
mass is - as above - the inverse of the second moment correlation function in the low 
temperature phase, 

m R = , (76) 

and hence vanishes at the critical point. 

One may keep in mind that we are not interested in the solution of this model, in a mass 
spectrum etc., but only in its behaviour near the critical point. It is the great advantage of 
universality that we are spared a detailed comparison between <I> 4 -theory and a statistical 
model of the binary fluid system. In order to make contact with experiment, one has to 
extract the analogue of the interface tension a in $ -theory. A suitable definition comes 
from considering tunneling in a finite volume. Consider a cylinder-type geometry, where 
the Euclidean space is a square of side-lengths L and periodic boundary conditions, while 
the Euclidean time X3 remains non-compact. The Hamiltonian of $ 4 -theory in 2 + 1 
dimensions is the generator of translations in the X3- (i.e. time-) direction. In the infinite 
volume limit L — > 00, the symmetry $ — > — $ of the Hamiltonian is spontaneously broken 
at low temperatures, since the field acquires a nonzero vacuum expectation value 

($) = ±v , (77) 

and two different but energetically degenerate ground states 

(+|$|+) = +v or (-|$|-) = -« (78) 
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exist. On the other hand spontaneous symmetry breaking does not occur in a finite 
volume: Rather, the degeneracy between the two ground states is lifted by tunneling 
effects. There exists a unique vacuum \S), which has energy E s = and is symmetric, 

(S|$|S)=0, (79) 

and another antisymmetric state \A) with an energy Eqa, which in the infinite volume 
limit approaches zero. This is of course exactly the result of the WKB approximation 
of the double well potential in ordinary quantum mechanics in the case that the wall 
between the wells is very high. Sidney Coleman's presentation of the double well in 



field theory is still unbeaten, and for details of the following presentation, the reader may 
consult his lectures. 

Symmetric and antisymmetric state are in the finite volume to lowest order given by 

\S) = ^(l+> + h» 

\A) = -^(|+)-|-». (80) 

The tunneling rate is given by the transfer matrix sandwiched between the states ( +| and 
|± ), T denoting here and in the following Euclidean time: 

-ht\ ■ \ _\(^ < o -e 0A t\ 



(+\e- HT \±) = -(l±e-^ 1 ) . (81) 

I will now substantiate that the energy of the lowest antisymmetric state vanishes expo- 
nentially with increasing L, and that it is related to the interface tension by (see Refs. |T5| ], 

B 0) 

E 0A = Ce~ aL2 . (82) 
The transition amplitude in the functional integral formulation is 

Urn (+|e-" T |-> = J V<§> e~ 5[ * ] (83) 

with the boundary conditions 




$(x) - U(h ^2 ■ (84) 



In the semiclassical (WKB) approximation the functional integral is expanded around the 
configurations of least action. The classical kink solution to the equations of motion, 



, ^^tanh^(x 3 - a) , (85) 
\ 9o 2 



has the smallest action of all configurations interpolating between +vq and —vq. Here, a 
is the free parameter which specifies the location of the kink on the 23-axis. The classical 
energy of the kink is related to the interface tension, 



S c = 2^L 2 = a L 2 , (86) 
9o 
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Figure 4: A typical, dominating 3-kink-contribution to the transition amplitude 
fl83|) of the matrix element ( +\e~ HT \— ) at a temperature just below the critical 
one. Shaded regions show the location of the kink, the ±-signs in between the 
adopted vacuum state. 



in this approximation since the reduced interface tension a of the statistical system is 
exactly the free energy of the kink per unit surface, i.e. its action per unit surface in the 
Euclidean field theory. 

As shown in Fig. f|, one has to take into account that the system can tunnel several 
times between the two states |+ ) and |— ), and in the WKB approximation, each of these 
tunneling processes will (when only tunneling takes long enough because of a sufficient 
barrier height) contribute to the functional integral a factor 

Ke~ Sc T . (87) 

The factor K can be calculated, but I shall not discuss it here. 

Therefore, having to cross an even (odd) number of interfaces in time development 
from a state | + ) (|— )) to | + ), one can estimate (|81~D by 



+ \e- HT \±) ~ Y K n e- nS *i- 

n even(odd) 



= exp (ife- 5c T) ± exp (-Ke~ Sc T) . (88) 

The factor T n /n\ comes from integration over all different locations of the kinks, i.e. taking 
into account double counting. Comparison with (|8lJ) indeed shows now (82) , 

E 0A = Ke~ Sc = Ce~ aL2 (89) 

in 3 dimensions, so that one knows now how to obtain the interface tension in the finite 
volume: It can be calculated from the energy splitting between the ground state and the 
first excited state. For L — > oo, E a — > and spontaneous symmetry breaking is possible, 
all tunneling processes being exponentially suppressed. In the following, we adopt (|89| ) 
as definition of the interface tension in field theory. 



3.3 Calculating R + in Field Theory 

This section outlines the work of references |18| and |]l9fl . In three-dimensional field theory, 
the coupling go has a positive mass dimension and - albeit the theory is hence super- 
renormalisable - infrared divergencies forbid one to construct the critical, i.e. massless, 
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theory perturbatively. One way out is to consider the theory in 4 — e dimensions and 
extrapolate to the case e = 1, but it was this technique which led to the unacceptable 
value of R + = 0.20 in an (9(e 2 )-calculation [13]. Also it turned out that convergence was 
very poor, the C(e)-calculation giving R + = 0.14. 

Another technique flll| |, ma de renormalisation group calculations in d = 3 possible 
by renormalised perturbation theory: Adopting a renormalisation scheme in which the 
renormalised coupling is denoted g R , the perturbative expansion goes in terms of the 
renormalised, dimensionless variable 



9r_ 

m R 



(90) 



Renormalised expansions in terms of ur do not show an infrared problem since even at 
the critical point, where m R — > 0, ur — > u* R remains finite. Thus, information about the 
critical theory can be obtained by perturbation theory at u* R , but ur is not a parameter 
but fixed by the fixed point value u* R , neither is it small. 
In our case we choose the renormalised coupling as 



9r 



m 



R 
?R 



(91) 



The fixed point value of u R is then 



u , 



15.1 ± 1.3 



(92) 



from various analytical results in the literature. 

Clearly, as one approaches the fixed point, i.e. the point of transition between mixed 
and separated phase, the interface tension vanishes and quantum fluctuations should 
become important. They are the large entropy fluctuations of statistical mechanics. The 
saddle point approximation of the functional integral takes into account the quadratic 
fluctuations about the classical kink solution (|85|) 



$(x) 



with 



M 



$ c (x) + 7j(x) 
S = S r +- I d 3 X 



s c + - 

2 



r/(x) Mi)(x) + 0{r] 3 



-ml cosh 2 ~^-{x3 



(93) 



(94) 



and corresponds to a one-loop perturbative calculation. 

The (lengthy) calculation is for the two-dimensional case outlined in Coleman's lectures 
i~4jl , and in the three-dimensional case it has been performed in Ref. [HJ. The result is 



that the fluctuations modify the energy of the antisymmetric state to 



OA 



'SI' 2 * 
2tt 



det'M 



detM 



-1/2 



where 



Ma 



(95) 
(96) 
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and all contributions from multi-kink configurations have again been taken into account. 
Because M has zero eigenvalues corresponding to translations of the kink in time (pa- 
rameter a), the determinant of M is calculated without the zero modes, as indicated by 
the prime over the determinant. The zero modes have to be treated separately by the 
method of collective coordinates J2T| and give rise to the factor S^ 2 oc L. 

The determinant can be calculated exactly for flat interfaces using heat kernel and (- 
function techniques |T£| . Three types of contributions of the final result are worth noting: 
First, it produces the counter-terms which are required to convert the unrenormalised 
parameters mo , go , Vq of (|85| ) into renormalised ones, giving secondly an additional factor 
1/L. Finally, it gives a substantial one-loop correction to the term proportional to L 2 in 



the exponent, and hence to the interface tension 

Omitting the details of the calculation of the determinant as well, one arrives by 
comparison with (|89| ) at an interface tension a, which has a negligible, exponentially 
small L-dependence and is in the infinite volume limit given by a Laurent series in ur, 
namely 



m 



a 



Ur 



47T 



15 



In 3 )+0(u 2 R ) 



(97) 



'39 

.32 16 

Comparison with the tree level result ( |S6"1) shows that the correction is not negligible. 
With (H/||/|9D and fl7|), one finally obtains 



) 2 = O"0 


(e - 


) 2 = 


«B f 


39 


15 


4tt V 


32 
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2 

U*R 

0.1024 ±0.0088. 



In 3) +0(u 



(91 



The one-loop contribution amounts to 22% of the leading term. 

As noted at the beginning, experimentalists have measured R + , i.e. above the critical 
temperature, while we calculated R below T c , but with the help of another universal 
quantity, the conversion factor 



^- = 1.96 ±0.03, 
to 



one finally obtains 



0.39 ±0.03 



(99) 



(100) 



in good agreement with the experimental value ([?[]). A more recent calculation of my 
Diplom- student Siepmann gave u* R = 14.73 ±0.14, leading to -R_ = 0. 1057 ± 0.0010. This 
number is in excellent agreement with recent Monte Carlo calculations by Hasenbusch, 
Pinn and coworkers 



23], which resulted in 



R. 



0.1056 ±0.0019. 



;ioi) 
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